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We study the pairing mechanism in iron pnictide superconductors based on the five-orbital 
Hubbard-Holstein model. Due to Fe-ion oscillations, the s-wave superconducting (SC) state without 
sign reversal (s++-wave state) is induced by orbital fluctuations by using realistic model parameters. 
The virtue of the present theory is that the famous empirical relation between Tc and the As-Fe-As 
bond angle is automatically explained, since the electron-phonon (e-ph) coupling that creates the 
orbital fluctuations is the strongest when the As4-tetrahedron is regular. The negative iron isotope 
effect is also reproduced. In addition, the magnitude of the SC gap on hole-pockets is predicted 
to be rather insensitive to the corresponding d-orbital {xz/yz- or 2;^-orbital), which is consistent 
with the recent bulk-sensitive angle-resolved photoemission spectroscopy (ARPES) measurement 
for (Ba,K)Fe2As2 and BaFe2(As,P)2. These obtained results indicate that the orbital-fluctuation- 
mediated s++-wave state is a plausible candidate for iron pnictides. 

PACS numbers: 74.70.Xa, 74.20.-z, 74.20.Rp 



I. INTRODUCTION 

The understanding of the pairing mechanism in iron 
pnictide superconductors [1] has been a significant open 
problem. By taking account of the Coulomb interaction 
at Fe-ions and the nesting of the Fermi surfaces (FSs), a 
fully-gapped sign-reversing s-wave state (s±-wave state) 
has been proposed based on the spin fluctuation theo- 
ries [2, 3]. Up to now, spin-fluctuation-mediated uncon- 
ventional superconductivity is believed to be realized in 
various metals, such as high-Tc cuprates [4-6], n-^BEDT- 
TTF)2X [7-9], and CeMIns (M=Co,Rh,Ir) [10]. To con- 
firm the spin- fluctuation scenario in iron pnictides, it is 
of significant importance to find evidences for the sign- 
reversal in the superconducting (SC) gap, and for the 
relationship between spin fluctuation strength and the 
SC transition temperature Tc. 

In principle, spin-fluctuation-mediated superconduc- 
tivity is fragile against nonmagnetic impurities or ran- 
domness, since the SC gap function has sign changes in- 
evitably. This is also true for iron pnictides, although 
the FSs are disconnected and the SC gap is fully-gapped. 
According to Ref. [11], decrease in Tc per pimp = 1/il^cm 
reaches ~ IK independently of the impruity potential 
strength. Contrary to this expectation, the SC state is 
very robust against various impurities [12] and heavy- 
particle irradiations [13, 14], although carrier number de- 
pendence may exist [15]. Moreover, the spin-fluctuation- 
mediated superconductors are expected to show a "reso- 
nance peak" in the neutron inelastic scattering as a re- 
flection of sign-change in the SC gap, as observed in high- 
Tc cuprates [16-18] and in CeMIns (M=Co,Rh,Ir) [19]. 
However, the observed "resonance-like" peak structure 
in iron pnictides [20-22] is reproduced theoretically by 
considering the strong correlation effect via quasiparticle 
damping, even in the conventional s-wave state without 
sign reversal (5++-wave state) [23]. 



In BaFe2(Asi_a^Paj)2, Tc increases as x decreases till 
the lattice structure transition occurs at x = 0.27, 
and Tc is positively correlated to the spin-fluctuation 
strength for x > 0.33 [24]. On the other hand, Tc 
in LaFeAsOi-^pFa; at x = 0.14 increases from 26 K to 
43 K by applying the pressure, whereas spin- fluctuation 
strength observed by 1/TiT measurement is almost un- 
changed [25]. Thus, the correlation between the spin- 
fluctuation strength and Tc seems to depend on com- 
pounds. 

Considering these difficulties in the 5±-wave scenario, 
we have proposed the orbital-fluctuation theory in Ref. 
[26], by taking account of the d-orbital degree of free- 
dom in iron pnictides. It was found that large or- 
bital fluctuations are induced by the electron-phonon (e- 
ph) interaction due to Fe-ion oscillations, although they 
are not induced by Coulomb interaction alone. Then, 
orbital-fluctuation-mediated s++-wave SC state is re- 
alized [26, 27] even if e-ph interaction is smaller than 
that estimated by the ffist principle study [28]. Exis- 
tence of large ferro-orbital fluctuations is suggested by 
prominent softening of shear modulus in Bal22 [29, 30]. 
Also, Raman spectroscopy [31], angle-resolved photoe- 
mission spectroscopy (ARPES) [32], and optical conduc- 
tivity measurement [33] highlight the importance of e-ph 
interaction. 

In this paper, we analyze the five-orbital Hubbard- 
Holstein (HH) model for iron pnictides in detail, by tak- 
ing account of all the matrix elements of the e-ph inter- 
action due to Fe-ion oscillations correctly. It is found 
that a small e-ph interaction (A < 0.15) can induce sub- 
stantial orbital fluctuations, utilizing all five d-orbitals 
on the FSs efficiently. Our main findings are as follows: 
(i) Empirical relation between Tc and the As-Fe-As bond 
angle (Lee plot) [34] is automatically explained, (ii) Ex- 
perimental negative isotope effect [35] is reproduced, (iii) 
Magnitude of the SC gap on the Z^-orbital hole-pocket 
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in (Ba,K)Fe2As2 and BaFe2(As,P)2 is comparable to the 
gap on other hole-pockets, which is consistent with exper- 
iments [36]. The range of model parameters for the 
wave SC state becomes wider in the presence of small 
amount of nonmagnetic impurities [11]. In addition, the 
predicted orbital fluctuations had been observed as the 
softening of the elastic constants C44 and C^; in Ref. [30]. 

Here, we consider the "orbital physics" in iron pnic- 
tides: It has been revealed that the ordered phase in 
mother compounds is not a simple spin- density- wave 
(SDW) state, but prominent orbital-density- wave coex- 
ists. In fact, recent bulk-sensitive ARPES in BaFe2As2 
below Tn had shown that the Fermi surface around F- 
point are mainly composed of xz-orbital, indicating the 
prominent nonequivalence between xz- and ?/z-orbitals at 
the Fermi level in mother compounds [37]. Also, appar- 
ent in-plane anisotropy of resistivity had been observed 
in detwinned Bal22 even above the structural transition 
temperature (Tg), suggesting the existence of nematic or- 
der as a pure electronic origin [38]. These experimental 
facts indicate the existence of orbital fluctuations even 
in (doped) superconducting compounds, and therefore 
orbital-fluctuation-mediated 5++-wave SC state is ex- 
pected to occur next to the orbital-density-wave state 
in iron pnictides. 



II. MODEL AND HAMILTONIAN 

We construct the flve orbital HH model for iron pnic- 
tides, by adding the electron-phonon (e-ph) interaction 
to the Hubbard model in Ref. [2] . The Hubbard model is 
comprised of the potential term e^, hopping term t^^, in- 
traorbital Coulomb t/, interorbital Coulomb U' ^ Hund's 
coupling J, and pair hopping J'] [2] 

i jji a ij jjiv a 

i ^ fj, fi>u 
-J ^ Si^ • Si,j + J' ^ ^Att^/u^^^^^^^^t j ' (1) 



is given by [26], 



where i,j denote the sites, and \i^v are the flve d-orbitals. 
We denote Z^, XZ, YZ, - , and XY orbitals 
as 1, 2, 3, 4, and 5 respectively: The X- and F-axes 
are parallel to the nearest Fe-As bonds, and Z-axis is 
perpendicular to the FeAs plane [2] . The XF-coordinate 
is given by —45° rotation of the x^y-coordinate spanned 
by the nearest Fe-Fe bonds around the z-axis. 

Now, we derive the e-ph interaction term due to the 
Einstein oscillation of Fe ions. The Coulomb potential 
for the d electron at r (with the origin at the center of 
the Fe ion) due to the surrounding As^~-ion tetrahedron 
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(r;«) = 3e2^{|r 



u-m 



\r-Ri 



^ ±A [2XZ • ux - TYZ • UY + {X^ - Y'^)uz] ,(2) 

where u is the displacement vector of the Fe ion, 
Rf is the location of surrounding As ions in Fig. 
1 (a); When AS4 tetrahedron is regular tetrahedron 
{a = 109.47°), a = ^2/3 and b = y/l/3. That is, 
R-^/RFe-As = (±a,0,6) and (0, ±a, -6) for Fe^^^, and 
R-/RFe-As = (±a,0, -6) and (0,zba,6) for Fe^^^, and 
A = 30e^/v/3i?^g_As- i^Fe-As is the Fe-As bond length. 
We neglect the As-ion oscillations since they do not in- 
duce substantial orbital fluctuations unless very large e- 
ph interaction is assumed. 

Nonzero matrix elements = 
given as 



XYZ ^ 



^34 



^22 

4 = 



'13 



^14 



^35 



^25 



zb2AaV7, 
zb2AaV7, 
#3 = ±2AaV7, 
= =b2AaV7-(l/\/3), 
= T2AaV7-(l/v^), 
= T2AaV7-(2/V3), 



(3) 



where a is the radius of d orbital. The obtained e-ph 
interaction does not couple to the charge density since 
v^j^ is trace- less. Thus, the Thomas- Fermi screening for 
the coeflicient A is absent. We stress that there are many 
nonzero (off-diagonal) elements in Eq. (3) due to the fact 
that As ions locate out of the Fe plane. Then, the e-ph 
interaction term is given by 



-ph 



EEEE 



(4) 



which represents the orbital exchange process induced 
by the Fe ion displacement. We will show later that sub- 
stantial orbital fluctuations involving all flve orbitals are 
developed because of many nonzero elements in Eq. (3). 

Next, we derive the phono n- mediated electron-electron 
interaction. The local phonon Green function is 



D{uji) 



(5) 



which is the Fourier transformation of {TrU^{r)u^{0)) 
(/i = X, y, Z). ujd is the phonon frequency, and uq = 
^/h/2MFeL0D is the uncertainty in position for Fe ions; 

uo = 0.044 A for ujb = 0.02 eV. cui = 2ttIT is the boson 
Matsubara frequency. 

For both Fe^^^ and Fe^^\ the phonon mediated inter- 
action is obtained as 

jJLUjJi'v' i (T(t' 
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In Eq. (6), nonzero V^j^^^fjyf are given as 

V24,24 = "^4,34 = V22,22 = Vas^as = — V22,33 = —d{^l)i 
^25,25 = ^35,35 = V24,35 = "^25, 34 = " ^(^z), 
Vl2,35 "1^13,25 "1^12,24 -1^13,34 - (1/ V3)^(a;^) , 

^12,12 = ^13,13 = -(l/3)^(a;0, 

"^^14,33 = -"^14,22 = -(2/^/?>)g{uJl), 

^14,14 = -(4/3)^(c^0,(7) 
where ^(cj/) = {2Ao? /lYDiuji). The following relations 

hold; Vim.l'm' = Vml^l'm' = Vim,m'U and Vim,Um' = 

Vi'm'^im- In our earlier investigation [26], we have es- 
timated ^(0) ^ 0.4eV for R^e-As ^ 2.4 A, a ^ 0.77 A, 
and cjd ^ O.OlSeV. We also obtain ^(0) ^ 0.33eV for 
ujd ~ 0.02eV. 



is the chemical potential, and is the kinetic term of 
Eq. (1). Then, the susceptibilities for spin and charge 
sectors in the RPA are given by [39] 



r{q) 



1 - r«x° (q) 

X° (q) 



where 



h=l2= h 

h = h ^ h 
J' h = h 7^ h 



k 

u 
k 
h 



(9) 
(10) 



(11) 




a=sin(a/2) 
b=cos(a/2) 



(b) 



X 



s{c) 



X 
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FIG. 1: (Color online) (a) As4 tetrahedron in iron pnictides 
shown along the z axis. Here, we put Rfs-As = 1- The As- 
Fe-As bond angle a is illustrated in Fig. 6 (a), (b) The 
diagrammatic expression for x**^^'*- 

In our earlier calculation [26] , only the first line of Eq. 
(7) was taken into account since the weight of orbit als 1 
and 5 on the FSs is small. In this paper, however, we will 
show later that orbital fluctuations increase substantially 
if all the interactions in Eq. (7) are taken into account 
correctly. This is one of the main message in the present 
work. 

Now, we perform the RPA for i^Hub + ^^.e- The ir- 
reducible susceptibility in the five orbital model is given 

by 



X?r,W (q) = E ^irn (k + q) Gl,, (k) , (8) 



where G^{k) = [icn + /i — H^]~^ is the d electron Green 
function in the orbital basis, q = {q^uoi)^ k = (/c, e^), and 
€n = (2n + l)7rT is the fermion Matsubara frequency, ja 



r%u;i) = -C-2V{uJi), 



U, ll=l2 = h 

-U' + 2J, h=h^h 

2U' - J, li=l2^ h 

J', li=l-i^ h 



h 
h 
h 
h 



(12) 



(13) 



Here, we neglect the ladder-diagram for phonon- 
mediated interaction because of the relation ujd <C 
^band[26, 27]. The Bethe-Salpeter equation for x^^^^ is 
given in Fig. 1 (b). 
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FIG. 2: (Color online) Obtained U-g{0) phase diagram for 
n = 6.1. Near the orbital-density- wave boundary, s++-wave 
SC state is realized by orbital fluctuations. 

Hereafter, we assume that J = J' and U = /7'+2 J, and 
fix the ratio J/U = 1/6. Figure 2 shows the U-g{0) phase 
diagram for n = 6.1 given by the RPA. ^^(c) is the spin 
(charge) Stoner factor, which is given by the maximum 
eigenvalue of r^^'^'^x^{q^O). The transition line for the 
spin (orbital) order is given by the condition as(^c) — 1- 
Note that this phase diagram is independent of ud since 
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as(c) is free from lud. For U = 1 eV, the critical value 
^cr(O) for Qfc = 1 is 0.22, which means that dimension- 
less coupling constant is Acr = gcr{0)N{0) ~ 0.15, where 
N{0) is the density of states per spin at the Fermi level. 
Hereafter, the unit of energy is eV. 

In our earlier work [26], we have shown that ^cr(O) ~ 
0.4 when only the first line of Eq. (7) is considered. How- 
ever, we stress that ^cr(O) in Fig. 2 in the present paper 
is almost halved. This result means that strong orbital 
fluctuations are induced by much smaller ^(0) by utiliz- 
ing all flve d-orbitals on the FSs efficiently. Therefore, we 
take all the interaction in Eq. (7) into account in later 
calculations. 




(c) 




(TT.O) 



FIG. 3: (Color online) Obtained (a) X24,24(q^^) and (b) 
X22,22(<7,0), for n = 6.1, [/ = 1, T = o'o2, and ^(0) = 
0.21. Note that X24,24(<7,0) = XmMq'^^) and X22,22(<?,0) = 
X33,33(^^0), where is given by the rotation of q by 7r/2. 
In (c), charge susceptibility 0) = ^i^^Xii,mm{Q,^) is 

shown. We use 2048 Matsubara frequencies. 

Figure 3 shows the obtained Xw mm'i^^^) 
{IV, mm') = (24,24) and (22,22), respectively. Used pa- 
rameters are n 6.1, /7 = 1, T = 0.02, and ^(0) = 0.21, 
which correspond to ac = 0.98. X24 24 (X34 34) 
X22 22 (Xss 33) most divergent channels for elec- 

tron doped case. First, we discuss X24 24 i^ ^^iS- ^ (^)* 



has the largest peak near Q = (7r,0), which comes from 
the nesting between FS3,4 and FS1,2 in Fig. 4 (c), and 
the multiple scattering by ^24,24- Its second largest peak 
near gr = (0, 0) originates from the forward scattering by 
^24,24 in the FS3 or FS4 that is composed of 2-4 orbitals. 
X24 24(^50) X24 24(^50) comparable, meaning 
that the ferro-orbital and antiferro-orbital orders are 
highly frustrated. Note that X24 24(^5 0) ^ X24 24(0, 0) if 
we consider only the flrst line of Eq. (7) as demonstrated 
in Ref.[26]. 

We also discuss X22 22 i^ ^^ig- ^ 0^)- ^^s the largest 
peak near Q = (7r,0), which arises due to nesting be- 
tween FS3,4 and FS1,2. X22 22 ^iso has a slightly lower 
peak around (0, 0) that originates from the forward scat- 
tering in each FS. We note that the large enhancement in 
X24,24 and X34,34 (X22,22 ^nd X33,33) is caused by in-plane 
(out-of-plane) Fe-ion oscillations. On the other hand, the 
total charge susceptibility x%Q^ 0) = Y.i,m Xii^mmiQ^ 0) is 
not enhanced as shown in Fig. 3 (c), because of the re- 
lation X22,22(<7^0) ^ -X22,33(<7^0) P^]. Therefore, the 
origin of the superconductivity in the present model is 
not charge fluctuations, but orbital fluctuations that can 
develop without cost of the Coulomb potential energy. 

Finally, we discuss the softening in the elastic con- 
stants due to orbital fluctuations. Recently, Yoshizawa 
et at. have observed large softening in C44, Cee, and 
Ce in Ba(Fe,Co)2As2 [30]. The corresponding strains 
for C744, ^66, and Ce are e^z, exy, and exx - ^yy, 
respectively. Using the point-charge model, one can 
verify that the strains e^y (/i, = X, F, Z) induce the 
quadrupole potential on each Fe ion; (\)^y ex iiv. The 
corresponding matrix elements = and 

^\(\>xx - (\>YY\m) = (l)f^~^^ are proportional to \v^\ 
and \vf^\, respectively, where vf^ are described in Eq. 
(3). In the linear response theory, the enhancement 
in C^^ (^E^) is proportional to the quadrupole susu- 
ceptibility E^z'w^lx?r,w(O'0)^w with E = XZ 
(S = - F^). As shown in Fig. 3, both X24,24 
and X22 22 evolve in the present study. Considering 
that (X \v^\ is flnite for /,m = 2,4, we find that 

the softening in C44 is induced by fluctuations in the 
(2, 4)-channel. In the same way, the softening in Ce 
is induced by fluctuations in (2,2)-, (3,3)-, and (2,3)- 
channels. These results are consistent with the reports 
in Ref. [30]. That is, theoretically predicted orbital fluc- 
tuations in Figs. 3 (a) and (b) had been conflrmed ex- 
perimentally. 



III. ELIASHBERG GAP EQUATION 

In this section, we analyze the following linearized 
Eliashberg equation using the RPA by taking account 
of both the spin and orbital fluctuations on the equal 
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footing [39]: 



k' ,mi 



(14) 



si^ )^jri4Jn3( ^ )) 



where, 



for the singlet states. A^; is the eigenvalue of the gap 
equation, which approaches unity as T ^ Tc. Hereafter, 
we use 64^ k meshes, and 1024 or 2048 Matsubara fre- 
quencies. We perform the calculation at relativity high 
temperatures (T > 0.02) since the number of meshes is 
not enough for T < 0.02. 
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FIG. 4: (Color online) (a) (b) Obtained SC gap functions 
for (a)[/ = 1.09 and {h)U = 1.10, respectively. We put 
g{0) = 0.21 (ttc = 0.98), T = 0.02, and ujd = 0.02. They 
are normalized as iV"^ 1]^,^^ \Aim{k)\'^ = 1. We use 2048 
Matsubara frequencies, (c) FSs in the unfolded Brillouin 
zone. FS1,2 (FS3,4) are composed of 2,3-orbitals (2,3,4- 
orbitals). FS5 is composed of 1-orbital in (Ba,K)Fe2As2 and 
BaFe2(As,P)2, which will be discussed in Sec. V A. Note that 
FSS moves to (0, 0) in the folded zone. 

Figure 4 (a) and (b) show the SC gap on the FSs in 
the band representation for U = 1.09 and U = 1.10, re- 
spectively. We put n = 6.1, T = 0.02, ujd = 0.02, and 



^(0) = 0.21, which correspond to ac = 0.98. The hor- 
izontal axis is the azimuth angle for k point with the 
origin at r(M) point for FS1,2 (FS4). For U = 1.09, 
the 5++ state is realized by orbital fluctuations [26]. On 
the other hand, the s± state is realized for U = 1.10 
since the spin fluctuations dominate the orbital fluctu- 
ations. In this case, the boundary of the 5++ s± 
phase transition is uniquely defined since the obtained 
gap functions are always full-gap. If one introduces low 
concentration of impurities, the 5++-wave state is real- 
ized even for U > 1.1 [26], and moreover, the transition 
becomes full-gap s++ ^nodal s full-gap s± [40]. For 
a quantitative study of the line nodes observed in several 
122 compounds [41], the 3- dimensionality of the FSs may 
be indispensable. 

Figure 5 shows the obtained A^; for (a) = and (b) 
= 1, respectively, at T = 0.02. In case (a), A^; exceeds 
unity for ujd = 0.02 when 1 — ac — 0.02. A^; increases as 
ujd increases. Considering the relation ujjj (x 1/a/M, this 
result means the positive isotope effect. In case (b), A^; 
is larger than unity for ujd = 0.02 if 1 — o^c — 0.01. We 
stress that A^; decreases diS ujjj increases, which means 
that the negative isotope effect is realized. 

Now, we discuss the origin of the negative isotope effect 
in case (b). In the BCS theory, Tc in a single band model 
is given by [42] 



1.13ujd exp 



A-/i* 



(16) 



where A = gN{0) is the dimensionless coupling constant, 
and /i* is given as 



1 + /iln(e/a;D)' 



(17) 



where /i = UN{0). ji* is called the Morel- Anderson 
pseudo-potential. In general, ja^ <^ ja since the limit of 
Coulomb interaction e is much larger than ujjj. Now, we 
derive the coefficient /3 = -dlnTc/dlnM (Tc ex M"^), 
where M is mass of Fe ion. By differentiating Eq. (16) 
by M using the relation ujd 1/>/M, we obtain 
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*2 



(A-M*r 



(18) 



The value of (3 decreases from 1/2 as /i* increases, and 
becomes negative when fi* is larger than A/2. This can 
be realized when U is relatively large. Therefore, the 
negative isotope effect {j3 < 0) in Fig. 5 (b) is caused 
by the reduction in /i*. In other words, negative isotope 
effect originates from the enhancement of the retardation 
effect. 

Note that Eq. (18) is valid only for one-band model. In 
order to obtain the correct coefficient /3 in iron pnictides, 
we have to analyze the five-orbital model. 

In case of Fig. 5 (b), A^; exceeds unity only when 
ac ~ 0.99. However, A^; can reach unity for smaller ac 
when the temperature is much lower than cujj = 0.02. 
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(a) 




c 



FIG. 5: (Color online) Obtained A^; as function of 1 — ac 
for (a)[/ = and (b)[/ = 1, respectively, at T = 0.02. The 
former (latter) corresponds to the positive (negative) isotope 
effect. We use 2048 Matsubara frequencies. Inset in (a): T- 
dependence of A^; for ac = 0.98. 



To explain this behavior, we analyze the following single 
band gap equation [43]: 

X^A = T^^^V{zen). (19) 

k,n ^ 

Here we consider the BCS approximation V{k^ icn) = 
—gO{ujD — \^n—7TT\). After carrying out the k summation, 
we obtain 

A^A = 2^7V(0)^^A, (20) 



where (2nc + l)7rT = ujd- When ud '> T (i.e. ric > 1), 
Eq. (20) is solved as [43]: 

A^; = gN{0) ln(1.13a;D/T). (21) 

Therefore, A^; diverges logarithmically at low tempera- 
tures. On the other hand, when ujd — Eq. (20) is 
solved by putting nc = as 

A^; = 2^7V(0). (22) 

In this case, A^; does not depend on T. Inset of Fig. 5 (a) 
shows the T-dependence of A^; for = and ^(0) = 0.23 
{ac = 0.98). For uod = 0.20, A^; increases at low tem- 
peratures, in accordance with Eq. (21). In contrast, 
A^; slightly decreases for cj^ = 0.02 at low temperature, 
which might be due to the smaller size of k- or cj-meshes. 
Unfortunately, we can not perform the calculation below 
lOOK since the numbers of k- and cj-meshes are not suffi- 
cient. However, even if a;/:) = 0.02, A^; will increase below 
T ~ 20K. Therefore, A^; is expected to exceed unity at 
low temperatures even if ac <C 0.99. 

Finally, we briefly discuss the effect of self-energy E, 
which has been dropped in the present study. The quasi- 
particle damping 7 (= imaginary part of S) reduces both 
as(^c) and Ae- However, the dimensionless coupling con- 
stant for ac = 0.98 in the FLEX approximation is only 
A = ^(O)Ar(O) - 0.2 [40]. Thus, the present orbital- 
fluctuation scenario is justified even if the self-energy 
correction is taken into account. In the FLEX approx- 
imation, ^cr(O) for orbital-density- wave state and Ucv 
for SDW state become infinity in 2-dimensional systems 
since Mermin- Wagner theorem is satisfied [44]. 

IV. CRYSTAL STRUCTURE AND Tc 

In previous sections, we assumed that the As4 tetra- 
hedron forms a regular tetrahedron. However, it is well- 
known that As-Fe-As bond angle a, which is shown in 
Fig. 6(a), closely relates on Tc experimentally [34]. Here, 
we extend the theory for general bond angle a. First, we 
derive the e-ph interaction for general a. In this case, 
the potential ^(r; u) in Eq. (2) is not changed except for 
A] A is changed as 

306^ / 3V^ . na a\ 
^ = sin^ - cos - . 23 

When the As4 tetrahedron forms the regular tetrahe- 
dron (a = Q/Q = 109.47°), the term in the brackets 
in Eq. (23) takes the maximum value; unity. Thus, 
the effective interaction for a is given by ^(0,a) = 
^/(0)[(3V3/2)sin^(a/2)cos(a/2)]2, where g({)) = g{0,ao) 
In Fig. 6(b), we show ^(0,a)/^(0) as a function of a. 
When a deviates from ao, A'^ decreases rapidly. 

Figure 7 shows the obtained X22 22(^50) 
X24 24(^5 0) ^ — 120°, which corresponds to LaFePO. 
We' put ^(0) = 0.21, = 1, and T = 0.02, which are 
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109.47 
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FIG. 6: (Color online) (a) The definitions of the As-Fe-As 
bond angle a. (b) g{0, a) /g{0) as a function of a. 




(TT,0) 



FIG. 7: (Color online) Obtained (a) X24 24(9,0) and (b) 
X22,22(9,0) for a = 120°, ^(0) = 0.21, [7=1, and T = 0.02. 
We use 1024 Matsubara frequencies. 
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FIG. 8: (Color online) a-dependence of A^; for (a)[/ = and 
^(0) = 0.23, and (b)[/ = 1 and ^(0) = 0.21, respectively. 
We put T = 0.02, and ujd — 0.02. In (a), s++-wave state is 
always realized. We use 1024 Matsubara frequencies. 



equivalent to those in Fig. 3. Here ac = 0.98 is satis- 
fied when As4 tetrahedron is regular. Compared to Fig. 3 
(a), magnitude of X24 24 ^^^S- ^ (^) is largely suppressed 
due to the reduction in the phonon-mediated interaction 
g{0,a) shown in Fig. 6 (b). 

Figure 8 shows the a-dependence of A^; for (a)/7 = 
and ^(0) = 0.23, and {h)U = 1 and ^(0) = 0.21, respec- 
tively. In both cases, ac = 0.98 is satisfied when As4 
tetrahedron is regular. In case (a), the 5++ state is al- 
ways realized. In case (b), the 5++ state is realized for 
100° < a < 120°, whereas s± is realized for a < 95° or 
125° < Of, because the orbital fluctuations become infe- 
rior to spin fluctuations. In both (a) and (b), A^; for s++ 
state rapidly decreases when bond angle a deviates from 



ao. This result is consistent with the well-known experi- 
mental relationship between bond angle a and Tc [34, 45], 
and supports the realization of the orbital-fluctuation- 
mediated 5++-wave state in iron pnictides. 

Finally, we note that the e-ph interaction due to As ion 
oscillation does not take the maximum value at a = ao* 
For example, e-ph interaction by Ai^-mode, which is 
given by the change in the parameter b in Fig. 1 (a), will 
monotonically increase as a decreases. Therefore, the ex- 
perimental relation between Tc and a strongly suggests 
that the iron pnictides are not conventional BCS super- 
conductors due to charge fluctuations by Ai^-mode, but 
are the orbital-fluctuation-mediated superconductors due 
to Fe ion oscillations. 
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V. DISCUSSIONS 

In previous sections, we have analyzed both 5++- and 
5±-wave states based on the five-orbital HH model. Here, 
we discuss both states in more detail, by making compari- 
son between theoretical results and experimental reports. 



A. SC gap in the Z^-orbital hole-pocket 

Recently, bulk-sensitive ARPES measurements had 
been performed in (Ba,K)Fe2As2 and BaFe2(As,P)2 [36]. 
The observation have revealed the Z^-orbital hole-pocket 
around (0, 0) at /c^ ~ tt in the folded Brillouin zone ((tt, tt) 
in the unfolded Brillouin zone; see Fig. 4 (c)), which was 
predicted by the first principle LDA study around kz = tt 
for larger bond angle a [46]. Moreover, it was found 
that the magnitude of the SC gap in the Z^-orbital hole- 
pocket is as large as that in other hole-pockets composed 
of XZ/YZ- (and - Y'^-) orbitals. However, the SC 
gap will strongly depend on the orbital nature of the FS 
parts in spin fluctuation mechanism, since the Z^-orbital 
does not participate in the nesting [36] . 

Therefore, this "orbital-independent SC gap in 
(Ba,K)Fe2As2 and BaFe2(As,P)2" is a very crucial test 
for theories to understand the pairing mechanism. Here, 
we shift the Z^-orbital level by +0.32eV in the present 
model to reproduce the Z^-orbital hole-pocket at /c^ ~ tt 
in the 3-dimensional model [46], and analyze the orbital 
dependence of the SC gap in detail. 

Figure 9 (b) shows the 5±-wave SC gap functions ob- 
tained for U = 1.0 and ^(0) = 0. The obtained parameter 
is Ae = 0.37. (Note that Ae for 5±-wave state decreases 
when Z^-orbital hole pocket appears [2].) As we can 
see, spin-fluctuation scenario predicts very small SC gap 
on the Z^-orbital FS, since the spin correlation between 
electrons in Z^-orbital is very small: In iron pnictides, 
spin fluctuations due to the nesting are mainly induced 
by XZ/yZ-orbitals via intraorbital Coulomb interaction 
U between opposite spins. However, spin correlation be- 
tween different orbitals is much weaker, since the Hund's 
coupling J is much smaller than U. For this reason, the 
SC gap in the Z^-orbital hole-pocket is very small. Sim- 
ilar "orbital dependent SC gap" is considered to be real- 
ized in p-wave superconductor Sr2Ru04 [47]. 

On the other hand, strong orbital correlation exists 
for all d-orbitals in the present orbital-fluctuation sce- 
nario, since the e-ph interaction due to Fe-ion oscillation, 
Eq. (3), possesses many nonzero interorbital matrix ele- 
ments. For this reason, the SC gap in the Z^-orbital hole- 
pocket can be large. Figure 9 (c) shows the 5++-wave 
SC gap functions obtained for U = and ^(0) = 0.20 
(ofc = 0.98), induced by orbital fluctuations. The ob- 
tained parameters are Ae = 0.99 and ac = 0.98. As 
expected, the SC gap on the Z^-orbital FS becomes com- 
parable with that on other FSs. 

Therefore, the small orbital-dependence in the SC gap 
in (Ba,K)Fe2As2 and BaFe2(As,P)2 [36] supports the 



present orbital fluctuation scenario. For a quantitative 
study of this issue, large 3-dimensionality of the FSs 
in (Ba,K)Fe2As2 and BaFe2(As,P)2 may be important. 
This is an important issue for our future investigation. 
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FIG. 9: (Color online) (a) Band structure (in the unfolded 
Brillouin zone) given by shifting the Z^-orbital level by 
+0.32eV. The hole-pocket around (7r,7r), FSS in Fig. 4(c), 
is composed of Z^-orbital. (b) s±-wave SC gap function for 
JJ — \X) and 5^(0) = 0. (c) s++-wave SC gap function for 
[/ = and ^(0) 0.19. 



B. bond angle, pnictogen height 

In Sec. IV, we have explained the experimental rela- 
tion between the As-Fe-As bond angle a and Ae in iron 
pnictides by assuming the s++-wave SC state mediated 
by the orbital fluctuations. On the other hand, Kuroki et 
al. [46] had studied the same issue based on the spin fluc- 
tuation theory: In 1111 compounds, the Z^-orbital hole 
pocket discussed in Sec. VA appears as the pnictogen 
hight z = i?Fe-As cos(a/2) decreases. Then, Ae for s±- 
wave state quickly decreases, since Z^-orbital does not 
contributes to the spin fluctuations. In this scenario, Ae 
monotonically decreases as z does; The decrease in Tc 
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for a < ao cannot be explained without assuming the 
accidental balance between pnictogen height effect and 
another opposite effect. 

In the present orbital fluctuation theory, Ae for 5++- 
wave state is rather insensitive to the appearance of the 
Z^-orbital hole-pocket since Z^-orbital also contribute to 
the formation of the orbital fluctuations as discussed in 
Sec. VA. Therefore, Ae or Tc will be mainly controlled 
by the bond angle a, which is consistent with the exper- 
imental report [34]. 

C. pressure effect on Tc 

We also discuss the pressure effect on Tc in iron pnic- 
tides. In LaFeAsOi-yF^;, Tc increases from 26 K to 43 K 
in overdoped sample {x = 0.14) by applying 3^4 GPa 
pressure [25]. In this case, both 1/TiT and normal-state 
resistivity p around Tc are rather insensitive to the pres- 
sure [25, 48]. Similarly, Tc increases drastically in FeSe 
under pressure [45, 49]. On the other hand, Tc quickly 
decreases for NdFeAsOi-^ and TbFeAsOi-^ under pres- 
sure, accompanying the decrease in the temperature de- 
pendence of p (i.e, the inelastic scattering) [50]. 

To understand the pressure effect on Tc, we would 
have to consider the change in the bandwidth Wband, 
in addition to the bond angle a. In NdFeAsOi-^y and 
TbFeAsOi_y, the reduction in the inelastic scattering in 
p under pressure suggests the suppression of spin/orbital 
fluctuations. This change is expected to originate from 
the increase in Wband, which drives the system toward 
weak-coupling regime. Then, reduction in spin/orbital 
fluctuations under pressure should make Tc lower. 

Next, we discuss the possible origin of the enhance- 
ment in Tc under pressure. In FeSe, Tc increases un- 
der pressure ~ 8GPa, whereas ao — a{> 0) slightly in- 
creases. Here, we emphasize that the change in the Fe-As 
bond length i?Fe-As would be the key parameter: Under 
pressure, the effective interaction due to e-ph coupling is 
given by ^(0) = (i^Fe-As/^Fe-As)^^°(0), where the suffix 
represents the quantity at ambient pressure. According 

to Ref. [49], i^^e-As/^Fe-As = 2.38A/2.30A at 8GPa, and 
thus ^(0) ^ (4/3)^^(0). This prominent enhancement in 
^(0) under pressure might be the origin of strong increase 
of Tc in FeSe under pressure. 

D. iron isotope effect 

In Sec. Ill, we have discussed the iron isotope effect 
based on the orbital fluctuation scenario. Experimen- 
tally, Liu et al. reported that the iron isotope coefficient 
j3 = -dlnTc/dlnM is - 0.35 for SmFeAsOo.ssFo.is and 
Bao.6Ko.4Fe2As2 [51]. However, Shirage et al. had re- 
cently reported the negative (or zero) iron isotope effect 
for the same compounds; /3 ~ —0.18 for (Ba,K)Fe2As2 
[35] and f3 - -0.02 for SmFeAsO^; [52]. The reason for 
the discrepancy is yet unclear. 



We have shown in Sec. Ill that the coefficient P 
changes from positive to negative as the Coulomb in- 
teraction increases. The first principle calculations had 
estimated that the ratio U /W^and for 122 compounds 
is larger than that for 1111 compounds [53]. Then, 
negative (zero) isotope effect reported for (Ba,K)Fe2As2 
(SmFeAsO^) does not contradict with the orbital fluctu- 
ation scenario. 

In the spin fluctuation scenario, P becomes positive 
(negative) when the intra-pocket phonon-mediated at- 
tractive interaction g{0;uji) is superior (inferior) to the 
inter-pocket one g{Q;cJi) [54, 55]. Therefore, sign change 
in /3 can be explained if the nature of e-ph interaction 
largely depends on the compounds. 



VI. SUMMARY 

In the present paper, we studied the five-orbital HH 
model for iron pnictides, and found that s++-wave SC 
state is induced by orbital fluctuations in the presence 
of small e-ph interaction (A < 0.15). Strong orbital 
fluctuations are induced by multiple scattering processes 
due to the e-ph interaction, involving all five c?-orbitals 
on the FSs. We stress that the second-order process 
alone, which is usually studied in conventional BCS 
analysis [28], can neither induce large orbital fluctua- 
tions nor high-Tc 5++-wave SC state. Roughly speak- 
ing, Tc in the orbital fluctuation theory would be given 
as Tc ~ cc;Dexp(-l/A*), where A* ^ A(l - ac)~^ is the 
enhanced coupling constant, and thus it is much larger 
than T^^^ - CJD exp(-l/A). 

The virtue of this theory is that we can also explain the 
following issues which remain unresolved within the spin 
fluctuation theory: (i) empirical relationship between Tc 
and the As-Fe-As bond angle (Lee plot), (ii) negative 
iron isotope effect in (Ba,K)Fe2As2, and (iii) orbital- 
independent SC gap in (Ba,K)Fe2As2 and BaFe2(As,P)2 
observed by bulk-sensitive ARPES measurement [36]. 
Recently, theoretically predicted orbital fluctuations in 
Figs. 3 (a) and (b) had been confirmed by the soft- 
ening of the elastic constants C44 and Ce [30]. These 
obtained results support the idea of the s++-wave state 
mediated by orbital fluctuations in iron pnictides, next 
to the orbital-ordered state in mother compounds. 

Finally, we list several significant future issues. The 
e-ph interaction due to corrective oscillations (e.g., half- 
breathing mode) might be important to increase the 5++- 
wave Tc. To make quantitative comparison between 5++- 
wave and s-t-wave states, study of the self-energy and 
vertex corrections for x^^^^ and the Eliashberg gap equa- 
tion is highly desired. The FLEX approximation would 
be useful for this purpose [40]. 
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Appendix A: Phase Diagram and Orbital 
Fluctuations for J/U > 1/6 



(a) 



orbital order JIU=\IA 




(TT,0) (TT,0) 



FIG. 10: (Golor online) (a) Obtained U-g{0) phase diagram 
for J/U = 1/4, 1/5 and 1/6. (b) Obtained X24,24(<7,0) for 
J/U = 1/6 and n = 6.1. (c) Obtained X22,22(c?,0) for J/U = 
1/6 and n = 6.1. 



In the present paper, we studied the multiorbital 
HH model for general parameters under the constraint 
J/U = 1/6. In fact, based on the first-principle calcu- 
lation, Miyake et al. had derived the averaged J/U as 
0.4eV/2.5eV=l/6.3 for LaFeAsO and 0.45eV/3eV=l/6.7 
for BaFe2As2, respectively [53]. Moreover, very small 
value of J/U < 0.1 is required to reproduce the small 
magnetic moment in the SDW state within the mean- 
field approximation [56]. However, J/U is expected to 
be larger for usual iron-compounds, and there is no con- 
sensus on the value of J/U in iron pnictides up to now. 
In case of J/U = 1/6, Coulomb interaction enhances the 
orbital fiuctuations when ^(0) is fixed, as shown in Fig. 
2. However, it is highly desired to study the orbital fiuc- 
tuation for general value of J/U. In this appendix, we 
discuss this issue using the RPA, and show that charac- 
teristic nature of orbital fiuctuations does not infiuenced 
by J/U. 



Figure 10 (a) represents the U-g{0) phase diagram 
given by the mean-field approximation for J/U = 1/4, 
1/5 and 1/6. For J/U = 1/4, the Coulomb interaction 
reduces the charge Stoner factor ac when ^(0) is fixed, 
indicating the suppression of orbital fiuctuations. The 
value of gcr{^) for o^c = 1 is 0.26 when U = 1.0; the 
corresponding dimensionless coupling A = gN{0) is only 
0.18. For J/U = 1/5, ac is almost independent of U. 
Therefore, the value of ^cr(O) increases with J/U: The 
reason is that U {U' = U — 2J) reduces (enhances) or- 
bital fiuctuations. As for the spin correlation, the value 
of UcT for as = 1 decreases with J//7, indicating the en- 
hancement of spin fiuctuations. 



Figure 10 (b) and (c) show the obtained X24 24(^50) 
and X22,22(<7.0) for J/U = 1/4, U = I, and ac = 0.98 
(^(0) ?^'0.26) at T = 0.02. Comparing with Figs. 3 (a) 
and (b), it is found that the orbital susceptibilities are 
almost independent of J/U < 1/4 for a fixed ac. Thus, 
the present orbital fiuctuation scenario for iron pnictides 
would be plausible for J/U < 1/4. 
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